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The linearized Navier-Stokes equations for a system of superposed immiscible compressible ideal 
fluids are analyzed. The results of the analysis reconcile the stabilizing and destabilizing effects of 
compressibility reported in the literature. It is shown that the growth rate n obtained for an inviscid, 
compressible flow in an infinite domain is bounded by the growth rates obtained for the corresponding 
incompressible flows with uniform and exponentially varying density. As the equilibrium pressure at 
the interface Poo increases (less compressible flow), n increases towards the uniform density result, 
while as the ratio of specific heats 7 increases (less compressible fluid), n decreases towards the 
exponentially varying density incompressible flow result. This remains valid in the presence of 
surface tension or for viscous fluids and the validity of the results is also discussed for finite size 
domains. The critical wavenumber imposed by the presence of surface tension is unaffected by 
compressibility. However, the results show that the surface tension modifies the sensitivity of the 
growth rate to a differential change in 7 for the lower and upper fluids. For the viscous case, the 
linearized equations are solved numerically for different values of Poa and 7. It is found that the 
largest differences compared with the incompressible cases are obtained at small Atwood numbers. 
The most unstable mode for the compressible case is also bounded by the most unstable modes 
corresponding to the two limiting incompressible cases. 

PACS numbers: 47.20.-k,47.20.Bp 

I. INTRODUCTION 

The Rayleigh- Taylor instability, which occurs due to the gravitational instability of a heavy fluid overlying a lighter 
fluid^*^ is of fundamental importance in a multitude of applications ranging from the turbulent mixing in Inertial 
Confinement Fusion 3,4 to astrophysical phenomena^*^ Small perturbations of the interface between the two fluids 
grow to large amplitudes. At early times, for small enough initial perturbations, the flow can be described by the 
linearized equations and the amplitude grows exponentially. Later, the interface evolves into bubbles of lighter fluid 
and spikes of heavier fluid penetrating the opposed fluid. If the initial interface is randomly perturbed then bubbles 
and spikes of different sizes are generated. Mathematically, the Rayleigh- Taylor instability is an ill-posed problem 
and the dependence on initial conditions is still an interesting question^ 

The theory for the linear stage for incompressible fluids agrees well with the experiments^^ In the absence of 
surface tension and viscosity, the growth rate increases indefinitely with the wavenumber. This trend is changed 
by the presence of viscosity, in which case the growth rate has a peak value and decreases towards zero for large 
wavenumbers. On the other hand, the presence of surface tension stabilizes perturbations with the wavenumber 
larger than a critical value. 

The role of compressibility on the development of the Rayleigh- Taylor instability between inviscid fluids has been 
studied by several au thorS)2*2ii iiLi2iiii4iiSii& however its effect compared with the incompressible case is still under 
debate. The earlier studies of the linear stage dealing with ideal fluids&iSi introduced simplifying assumptions and 
are strictly valid only when 7 = 1. Bernstein and BooltiS and Turner— removed these assumptions and studied the 
effects of compressibility as a function of 7. They found that these effects are more important at small wavenumbers 
and the rate of growth increases as 7 decreases. They concluded that compressibility has a destabilizing effect. The 
same conclusion is obtained for a multilayer system by Yang and Zhangi^ by comparing the compressible growth 
rate with that corresponding to an incompressible system obtained as 7 — > 00. The increase in the growth rate as 7 
decreases can also be explained using the energy principle, as a special case of the comparison theorem in the calculus 
of variations— 

On the other hand, Sharp 7 finds a stabilizing effect of compressibility. Moreover, numerical results for late time 
growth seem to indicate that a (the constant of proportionality in the quadratic law for the rate of growth) increases 
with the speed of sound^ so that compressibility would have a stabilizing effect. Baker— found both stabilizing 
and destabilizing effects of compressibility in the linear regime, however his results were based on previously derived 



2 



formulas using different assumptions than ideal gas. The role of compressibility on the instability growth is thus not 
yet settled. Moreover, to the best of our knowledge no study of the effects of compressibility in the linear regime 
for viscous fluids has been performed. Additionally, there is no systematic study of the effects of surface tension and 
finite size domain in the compressible case. 

In this paper we resolve the apparent contradiction between the stabilizing and destabilizing character of com- 
pressibility for ideal fluids. We show that compressibility can be characterized by two parameters, 7 and the speed of 
sound, with opposing influence on the instability growth. Moreover, as 7 or the speed of sound (varied by changing the 
equilibrium pressure at the interface while keeping the interface equilibrium density constant) increase to 00, the 
limiting incompressible flows (and the bases of comparison) are different. The compressible growth rate is bounded by 
the growth rates obtained for these two incompressible flows which have exponentially varying and constant density, 
respectively. 

A physical explanation for the decrease in the growth rate as p^, decreases can be formulated based on the influence 
of poo on the local Atwood number. As the interface develops, the heavier fluid reaches regions of larger and larger 
densities of the lighter fluid while the lighter fluid reaches regions of smaller densities of the heavier fluid. The local 
Atwood number in these regions away from the initial position of the interface depends on poc , since the equilibrium 
density profile depends on p^. Using the equilibrium density provided below, it can be shown that the local Atwood 
number is lower for the points on the interface above the initial position, while it is higher below, compared to a 
system at higher p^ (less compressible). However, the decrease in p^ leads to a larger change in the local Atwood 
number above the initial position of the interface. The overall effect would be a decrease in the average Atwood 
number, thus offering an intuitive argument for the decrease of the growth rate with decreasing p^ . Nevertheless this 
argument might break for small domain sizes and the validity of the results for finite size domains is discussed. The 
influence of surface tension and viscosity on the growth rate are also considered. 

II. LINEARIZED EQUATIONS 

The case of two superposed compressible ideal fluids separated by an interface at X\ = is considered. The fluids 
are subject to a constant gravitational acceleration g — (—5,0,0). For each fluid, the motion is governed by the 
continuity, momentum transport, and energy transport equations 

dp , Opu k = Q (l) 



dt dxk 
dpui dpUiUk dp dr., 



dt dxk dxi dxi 



pgSn (2) 
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where p is the density, Ui the velocity in n direction, p the pressure, e the specific internal energy, and T the 
temperature. The viscous stress is assumed Newtonian, = p ^ff 1 + §~ — f j ^jjjj > an d t ne kinematic viscosity, 
p, and thermal conduction coefficient, A, are considered constant. The equations 11131 should be supplemented with 
equations of state for the pressure and internal energy. For each fluid, the specific heats are assumed constant, 
p = RpT and e = c v T, where R is the gas constant and c v — R/(j — 1) is the specific energy at constant volume. 
With these assumptions e = p ^_^ an d the energy equation becomes 

dp duk dp t ! 1 x _ duj t fl 8 ( A dT 



dt ^ dxk k dxk ^ ^ J '* dxk ^ ^ dxk \ dxk ) ^ ^ 

A. "0"-order equations 

The two fluids are assumed initially at rest and the primary variables are written as small perturbations around 
the equilibrium ("0") state, with u = 0. Then the governing equations reduce to 
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dp 
dt 
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dpo dpo dp 

dx~, = - pa9 > 9^^ = ° (6) 
d Po , n d ( dT \ 



9i dx\ \ dx\ 

Moreover, the "0" t,l -order variables are assumed to be in steady state, so that ^f- = and Tq and po are continuous 
across the interface. For infinite domain or finite size domain in xi direction with adiabatic walls, the energy equation 
yields Tq — constant. Consequently, the equilibrium state is given by 

Pom = -5-^r eic P(—5-^- x V ( 8 ) 

P0 m = Poo6Xp(- 9 Xl) (9) 

To = constant, (10) 

where p^ is the unperturbed pressure at the interface (izq = 0) and m = 1, 2 denotes the material 1 or 2, with material 
2 chosen to be above material 1. 

B. First-order equations 

The interface between the two fluids is perturbed with an X2 and X3 dependent perturbation. The location of the 
interface is described by the function x s (xi, X2, X3, t), with = u±. Moreover, a surface tension is added at the 
interface between the two fluids. The first order linearized equations are obtained as 

^+p A + u 1 Dp = (11) 
at 

du\ d ( dui \ 2 

po ~dt = - Dp - p9+p dx-{d^ + Du >)~r DA 



d 2 x s 



dxidxo dx^dxr- 



S(xx-x s ) (12) 



du 2 dp d ( du-i duj \ 2 dA , . 



dt 8x2 dxj \ dxj 8x2 J 3 8x2 

^ du 3 = dp | 8 ( 8u 3 | duj \ 2 OA 
dt 8x3 dxj \ dxj 8x3 J 3 8x3 

dp 8 f 8T\ 

^ = - 7P0 A- Ul Dp + ( 7 -l)—[X— J, (15) 

where A = is the dilatation and D denotes Eg uat ions 1111151 should be supplemented by boundary conditions 
and jump conditions across the interface. These conditions will be discussed at length in the next sections. 

By examining equations II HI 51 it should be noted that in the absence of heat diffusion, the limit of incompressible 
flow (A = 0) can be obtained either by letting 7 — > 00 (as considered by the previous authors) or by letting — > 00 
(note that since the equilibrium density at the interface should not be affected by the change in pressure, the equation 
of state implies Tq — > 00 in this case). In the latter case the equilibrium density becomes po = constant, while in 
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the former the exponential variation is still allowed. Therefore the incompressible limit (and the base of comparison 
for the rate of growth) is different in the two cases. In the next sections 7 and are considered as independent 
parameters, both affecting the compressibility of the flow. The inviscid, infinite domain and no surface tension case 
will be examined first. Then the influence of finite size domain, surface tension and viscosity on the results obtained 
will be investigated. The non-zero heat diffusivity case was also considered, but the results were very close to the 
non-diffusive case and are not presented here. 



III. INVISCID CASE 



In general, the linearized equations do not admit analytical solutions. However, it is possible to obtain an analytical 
solution in the absence of viscosity and heat diffusion. Following the usual approach (e.g. Ref>£), we seek solutions 
whose dependence on X2, x 3 , and time have the form 

exp(i(k2X2 + k 3 X3) + nt) (16) 
where k2, k 3 and n are constants. For solutions having this dependence, equations II 11151 with p, = 0, A = 0, become 



np 



-p A - uiDpo 

-2r 



Pquui = —Dp — gp — kT s 5(xi—x s ) — 

n 

panu 2 = -ik 2 p 
p nu 3 = -ik 3 p 

np = --fpoA + Uipog, 



(17) 
(18) 

(19) 
(20) 
(21) 



where A = Du± + i(k2U2 + k 3 u 3 ) and k 2 = k\ + k 2 . After eliminating p, A, U2, and u 3 from the above equations, an 
equation for u\ is obtained as 
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= 0, 



(22) 



where c = y 7^ is the speed of sound. The jump condition at the interface can be obtained by integrating equation 
1221 over an infinitesimal element of x\ which includes the interface 



A) 



9/c 2 



k 2 + n 2 /c 2 



k 2 + n 2 /c 2 



fc2 t g x n 

— l s u s - u s —dpo = 0, 



(23) 



where 6f = /+ — /_, with /+ = f(x s + 0), /_ = f(x s — 0), is the jump of a quantity f across the interface. The 
subscript s denotes the value which a quantity, continuous at the interface, takes at x± = x s . 
On each side of the interface, c 2 and — are constant, so that equation 1221 becomes 



D 2 Ul 



-D Ul 



n 



( 7m - l)g 2 k 2 



ui = 0, 



(24) 



with the solution of the form u\ m = A m exp(Ai m xi) + B m exp(A2 m xi), where 




x h2 m = 7r^±kJl + - J - T + —^- 2 + -r 2 -r- (25) 



Formula [23 for A% 2 m wa s also obtained by Amala^ The coefficients A m and B m can be determined to a multiplying 
constant from the conditions that u\ vanishes at the rigid boundaries located at x\ = —l\ and x\ — I2 and that it is 
continuous over the interface. After replacing u\ in the jump condition 1231 a dispersion relation can be obtained as 
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n 2 = F t F 2 [a 2l2 { ll +n 2 Ma 1 ) - (72 + n 2 Ma 2 ) - T s (jx + n 2 Ma x ) (72 + n 2 Ma 2 )] 
^2 



/ [ai7i(72 + n Ma 2 )[Xi 1 exp(Ai 1 Li) - A 2l cxp(A 2l Li)]F 2 
- ^272(71 +n 2 Mai)[X l2 exp(-Ai 2 L 2 ) - A 22 exp(-A 22 L 2 )]Fi] , (26) 

where the nondimensional quantities are defined as n 2 = M = g( - < ^ +p2 - > , T 8 — g (J°'l P2 - ) 1 = kl m , Ai >2m = A \ 2 '" = 
SmM±,/l + ^n 2 a m A/ + + F x = expCXi.LiJ-exp^Li), and F 2 = exp(-A l2 i 2 )-cxp(-A 22 L 2 ). 



7m m 7m H 4 

The densities p\ = and p 2 = are the values of po on the two sides of the interface and a m — p m J (pi + p 2 ). 

Relation 1261 represents a generalization of the dispersion relations obtained by previous authors. It includes the 
effect of a finite size domain and surface tension, and makes no isentropic assumption. For L±,L 2 — > 00 the dispersion 
relation becomes 



_ 2 0272(71 + n Mai) — 0171(72 + n Ma 2 ) - Tjji + n Mai)h 2 + n Ma 2 ) . . 

n = — = — = . (Zi) 

0171(72 + n Ma 2 )X ll - o 2 7 2 (7i + n Mai)X 22 

As explained above, the incompressible flow limit (A = 0) can be obtained either by letting p^ — > 00 or 7 — > 00. 
As Poo — > 00, the flow approaches incompressible flow with constant density, for which the nondimensional rate of 
growth is given by 



—2 o 2 - 01 - T s 

n i = . 77 . : , - 7 , (28) 

Oi coth(_Li) + a 2 coth(L 2 ) 
2 

which for infinite domain becomes n 1 i = a 2 — ct\ — T s £ On the other hand, as 7 — > 00, the equilibrium flow has still 
exponential varying density, with the dispersion relation 



F X F 2 [a 2 - ai -T^/ a 1 [A J i 1 exp(A 4 i 1 ii) - A 4 2l exp(A l 2l ii)]F 2 l 
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where A l i, 2 „ 
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and arc defined using A*i, 2m . The parameter M is related to the 



exponent in the formula for the unperturbed density po m — p m exp(—Ma m kxi). For infinite domain and no surface 
tension equation 1291 reduces to the formula derived by Bernstein and Bool*^ 



n"1 = -Ma\Oi 2 + ^M 2 a\a\ + (a 2 - oti) 2 (30) 

2 2 2 2 

and it is easy to show that n" i < n' i for M > and n" \ — > n' \ as M — * 0. 

In the dispersion formula for the compressible case ( equation 126(1 the decrease in p^ is equivalent to an increase in 
M at 7 constant. Thus, for the compressible case, M represents a measure of the compressibility effects on the rate of 
growth. On the other hand, M is proportional to the ratio between the wavelength of the initial perturbation and the 
density exponential change length scale. For small values of M, this lengthscale is much larger than the wavelength 
of the initial perturbation, and the rate of growth approaches the incompressible, constant density result. 

For infinite domain, an approximate relation for n 2 , valid to order O(M) for small values of M is 
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(31) 



For 71 



72 



1 the relation derived by Plesset and Prosperettii* is recovered. In general, for 71,72 > 1, it can be 



shown that relation 1311 implies n 2 < n' \ for any combination of parameters, with the constraint p 2 > p\. Moreover, 
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Fig. 1 shows that n 2 decreases as Poo decreases. It should be noted that the same effects of decreasing p^ on the 
nondimensional growth rate can be obtained either by increasing g or decreasing k, with all other parameters kept 
constant. In other words, compressibility effects are more important at larger values of g and lower wavenumbers, 
which is supported by the numerical solutions of the dispersion relation presented in Fig. 1. 
For large M and no surface tension, it follows from the dispersion relation that 



- 2 « I — i (Q r Qi)2 (32) 

axa 2 (ax 2^1 + « 2 2i_I + J a* + a* + 2 ^%-*~>> +2 ai a 2 ) 

A further reduction in the rate of growth can be obtained by increasing the adiabatic exponents. At the limit, when 
7i:72 ~>- oo, the nondimensional rate of growth becomes 



_ 2 1 (a 2 - ai)" 



rT~— v - ' , (33) 

which is the same as that obtained from formula 1201 for large M by Bernstein and Book 1 — and Turner^ In general, 
for finite values of M it can be showni 2 *^ that increasing the ratio of specific heats leads to a decrease in the rate 
of growth, also supported by Fig. 1. However, the rate of growth has different sensitivities to the change of 71 and 
72. Thus, as Fig. 1 shows, the change in the ratio of specific heats of the lower fluid leads to a larger change of n. 
Therefore, the rate of growth is more sensitive to the change in compressibility of the lower fluid. However, as Fig. 1 
indicates, these results are also sensitive to the value of the Atwood number. For large values of the Atwood number 
the results obtained for the compressible cases show little sensitivity to changes in the ratios of specific heats and the 
rate of growth obtained for the compressible case is close to the incompressible variable density result. Nevertheless, 
at small Atwood numbers and for large M, the values of 71 and 72 become important in determining the rate of 
growth. In addition, the relative difference between the growth rates obtained for the two limiting incompressible 
flows increases as the Atwood number decreases, so that compressibility effects are larger at small Atwood numbers. 

In conclusion, the instability growth rate for a compressible flow (n) in the inviscid, infinite domain and no surface 
tension case is bounded by the growth rates of the corresponding incompressible flows obtained for uniform (nj) and 
exponentially varying density (n"), so that n'( < n < n^. As p^ increases (so that the flow becomes less compressible), 
n increases towards n' i7 while as 7 increases (so that the fluid becomes less compressible), n decreases towards n". 



A. Influence of Finite Size Domain 



For the case in which the domain is bounded by rigid boundaries located at x = —l\ and x — 1%, the growth 
rates obtained for constant density incompressible flow and incompressible flow with exponentially varying density 
are given by equations |^H1 and respectively. Figure 2(a) presents the growth rate as a function of the nondimen- 
sional parameter M for different domain sizes. For domain sizes not very small compared to the wavelength of the 
perturbation, the nondimensional growth rate is still bounded by and n" and the rate of growth decreases as the 
domain size decreases. However, similar to the incompressible case, the decrease in n is less significant when L\ is 
decreased. Therefore, for L\ < L2 the growth rate varies more when p^ is changed, so it is more sensitive to the 
change in compressibility. 

For the extreme case when L2 << 1 (domain size small compared to the wavelength of the initial perturbation) 
and 71 « 1 it is possible, as Fig. 2(b) shows, that the compressible growth rate becomes larger than the constant 
density incompressible growth rate for M smaller than a critical value. Numerical solutions of the dispersion relation 

123 for a large range of parameters indicate that the curve n 2 can intersect only once the line n 2 = n' i . Therefore, an 
analytical condition for the existence of the overshoot can be found by letting M — > (in which case the dispersion 

relation simplifies considerably) and imposing n 2 > n' { . After some algebra one obtains 



ot\ coth LI + a.2 coth i 2 



coth L\ + coth L2 
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(34) 
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where $ m = L m coth 2 L m — coth L m — L m varies from to —1 as L m increases from to oo. Consistent with the 
numerical results, condition 1341 can be fulfilled only for small values of the domain size of the upper fluid and ratio of 
specific heats close to 1 for the lower fluid. 

B. Influence of Surface Tension 

The presence of surface tension tends to inhibit the growth of the instability. Moreover, for the incompressible case 
there is a critical wavenumber k c = [(p2 — pijg/T] 1 ^ 2 ^ so that the arrangement is stable for k > k c . By imposing n = 
in the dispersion relation 1261 it can be seen that the critical wavenumber remains the same as in the incompressible 
case. However, for T s ^ 0, the wavenumber appears as an explicit parameter in the dispersion relation 1261 so the 
variations of poo and k are no longer equivalent. Nevertheless, at each k, the nondimensional compressible rate of 

2 2 

growth obtained for infinite domain is still bounded by n' i and n" \ , as Fig. 3 shows. However, the lower limit is 
approached differently as 71 or 72 increase to 00. Thus, the variation in the compressibility of the lower fluid is more 
important at lower wavenumbers, while the variation in the compressibility of the upper fluid is more important at 
higher wavenumbers. 

IV. EFFECT OF VISCOSITY 

Consider the case of two viscous fluids, bounded by two rigid surfaces at x = —l\ and x = I2. Following the previous 
procedure, the linearized equations can be reduced to a single fourth order ordinary differential equation in u\ of the 
form 

A 4 D i u 1 + A 3 D 3 Ul + A 2 D 2 Ul + A x Du x + A ui = 0, (35) 

where the coefficients Ai are given in the Appendix for the compressible case and the two incompressible limiting cases. 
However, only for the uniform density incompressible case this equation has constant coefficients and an easily derived 
analytical solution. The boundary conditions for equation 1351 are m — at x — —l\ and x — I2, Ui and tangential 
viscous stresses continuous at the interface, and a jump condition found from the integration of the governing equation 
over the interface. For the compressible case, the condition that tangential velocities vanish at the rigid boundary 
can be written as A — Du\ = at x — —l\ and x = I2, while the continuity of U2 and 113 at the interface leads to the 
continuity of A — Dui. The divergence of the velocity fluctuations is given by 

A = B 3 D 3 Ul + B 2 D 2 Ul + B 1 Du 1 + B ui, (36) 

with the coefficients Bi given in the Appendix. The continuity of the tangential viscous stress over the interface can 
be written as 

S(fi[DA - D 2 Ul - k 2 Ul ]) = 0. (37) 

An expression for DA in terms of the derivatives of u\ is provided in the Appendix. The jump condition at the 
interface can be found by eliminating p, u 2 , and 113 from the momentum equation and integrating the resulting 
equation over an infinitesimal element of x\ which includes the interface 

6[(- P + ^D 2 )(A-D Ul )] + -5(nD Ul ) = ~[g( P 2 - Pi) - k 2 T s ]u s 
n n n z 

2k 2 

(rt!-A*i)(A-D«i)a (38) 

n 

For A = condition |2H1 reduces to the condition derived for the incompressible case in Refi^ In the Appendix an 
expression for D 2 A is provided. Since u\ can be found only to a multiplying constant, the boundary conditions are 
supplemented with the specification of u\ or one of its derivatives at one point inside the domain. Then equation 1351 
together with the boundary conditions described above form a closed set of equations from which u\ on each side of 
the interface and the rate of growth n can be determined. For all cases considered l\ and I2 are large compared to 
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the wavelength of the initial perturbation so the configuration is close to the infinite domain case. Equation 1351 was 
integrated on each side of the domain using a fourth order Runge-Kutta scheme. In order to determine n and U\ from 
the matching conditions at the interface, a multidimensional secant method (Broyden's method) was employed. 

Figure 4 presents numerical solutions of the viscous linearized equations for different Atwood numbers. Consistent 
with the previous results, the compressible rate of growth is bounded by the incompressible rates of growth obtained 
for uniform density and exponentially varying density cases. Moreover, its behavior is similar to the well known 
incompressible constant density result. It has a peak at some critical wavenumber and decreases towards zero as the 
wavenumber becomes large. However, the location of the peak is different compared to the incompressible case, with 
the highest difference at small Atwood numbers. Here we should note a qualitative difference between the results 
obtained for the constant density incompressible case and those obtained for the compressible and variable density 
incompressible cases. At small Atwood numbers, the critical wavenumber decreases with Atwood number for the 
constant density incompressible case, while it increases for the other two cases. Moreover, as shown in the previous 
sections, the largest differences in the inviscid rate of growth compared to the constant density incompressible case 
are obtained at small wavenumbers and small Atwood numbers. It is expected then that, for the viscous case at small 
Atwood numbers, the relative difference between the rate of growth obtained for the compressible and incompressible 
cases considered will be largest, also confirmed by Figs. 4(a) and 4(b). Although this difference should become very 
small as M approaches zero, Fig. 4(b) shows that, for small Atwood numbers, it persists at smaller values of M. 

V. CONCLUSIONS 

The effects of compressibility on the growth rate of Rayleigh- Taylor instability between two immiscible ideal fluids 
are examined in the linear regime. The results distinguish between the stabilizing and destabilizing character of 
compressibility. For infinite domains, the growth rate n obtained for the compressible case is bounded by the growth 
rates obtained for the corresponding incompressible flows with constant and exponentially varying density, and this 
result is not affected by the presence of surface tension or viscosity. For ideal gases with zero heat diffusivity, the 
limiting incompressible flow (defined by dui/dxi — 0) can be attained either by increasing the ratio of specific heats, 
7, or the speed of sound (varied by changing the equilibrium pressure at the interface at constant equilibrium density 
at the interface). The equilibrium density distribution for the limiting incompressible flow is different in the two cases. 
Moreover, the two parameters have opposing influence on the rate of growth. As the speed of sound is increased, the 
rate of growth increases towards the value obtained for the corresponding constant density incompressible flow, while 
as 7 increases, n decreases towards the value obtained for the corresponding incompressible flow with exponentially 
varying density. The presence of heat diffusion was also considered, but the results were very close to those obtained 
for the nondiffusive case and were not presented here. 

The equilibrium density for a compressible flow varies exponentially with x\ and depends on p^. Therefore, the 
local Atwood number changes as the interface moves away from the original position. Compared to a flow with a 
higher value of p^ (less compressible) , the local Atwood number decreases for the points on the interface situated 
above the initial position, while it increases for the points on the interface situated below the initial position. However, 
the change in the local Atwood number is larger above the initial position of the interface, so that the overall effect 
would be a decrease of the average local Atwood number. This offers an intuitive argument for the decrease of the 
growth rate as Poc decreases. Moreover, this argument suggests that the bubble velocity decreases, while the spike 
velocity increases for more compressible flows. On the other hand, as 7 decreases the fluids are more compressible, 
however the equilibrium density and pressure do not change. Therefore, as the heavier fluid moves towards regions 
of higher pressures, its volume decreases and the volume change is larger for more compressible fluids, so that the 
spike velocity decreases. Similarly, for more compressible fluids the bubble velocity increases. If the two fluids have 
different values for 7, it is shown that the growth rate is more sensitive to the change in the ratio of the specific heats 
of the lower fluid. However, at large Atwood numbers the rate of growth is little influenced by the values of 7I and 
72 and poo becomes the main compressibility parameter. In addition, it is shown that compressibility effects are more 
important at small Atwood numbers. 

For domains bounded by rigid surfaces, the compressible growth rate is still bounded by the two incompressible 
growth rates described above, except for the extreme case when the domain size of the upper fluid is small compared 
to the wavelength of the initial perturbation and 7 « 1 for the lower fluid. In this case, the compressible growth 
rate can become larger than the growth rate obtained for the corresponding constant density incompressible flow for 
values of the compressibility parameter M = g ^ 1 p +P2 ' ) smaller than a critical value. An analytical condition for the 
existence of this overshoot is provided. In general, the results show that the compressible growth rate varies more 
when the rigid boundary of the lower fluid is closer to the interface than the rigid boundary of the upper fluid, so 
that it is more sensitive to the change in compressibility. 

The presence of surface tension tends to inhibit the growth rate of the instability and for the incompressible case 
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there is a critical wavenumber above which the configuration becomes stable. It is shown that the value of this 
critical wavenumber is not affected by compressibility. For wavenumbers below this critical value the general result 
presented above remains valid. However, the presence of surface tension modifies the sensitivity of the growth rate to 
a differential change in the value of 7 for the two fluids. At smaller wavenumbers, the change in 7 for the lower fluid 
is more important for the variation of n, while the opposite holds true at higher wavenumbers. 

Numerical solutions of the linearized equations show that for viscous compressible fluids, the growth rate behaves 
in a manner analogous to the incompressible growth rate. It has a most unstable wavenumber and decreases towards 
zero at larger wavenumbers. Moreover, both the growth rate and the most unstable mode are bounded by the 
values obtained for the corresponding constant and variable density incompressible flows. For the constant density 
incompressible flow it is known that the most unstable mode moves to small wavenumbers as the Atwood number is 
decreased. The inviscid results presented in this paper show that the effects of compressibility are more important 
at small wavenumbers and small Atwood numbers. Consistent with these results, it is found that for viscous fluids 
compressibility becomes more important at small Atwood numbers. For small enough Atwood numbers, the difference 
between the compressible and incompressible growth rates will remain sizable at larger values of the equilibrium 
pressure. 

An interesting question raised by the results presented in this paper is if they remain valid in the nonlinear regime for 
single and/or multimode initial perturbation. Our preliminary numerical results seem to indicate a similar influence 
of Poo and 7 on the growth rate (and on the spike and bubble velocities) to that found in the linear regime. Moreover, 
even for large values of the equilibrium pressure so that the early time results are close to the incompressible flow 
results, the late time spike and bubble velocities become different than in the incompressible case. Another interesting 
question is about the range of the amplitudes of the perturbation for which the growth rate agrees with the linear 
theory prediction. Again, our preliminary numerical results seem to indicate that the range of validity of the linear 
assumption remains approximately the same as in the incompressible case. These results will be published elsewhere. 

This study was concerned with the effects of compressibility on the instability growth between immiscible fluids 
with uniform equilibrium temperature. It does not cover many of the configurations of practical interest, for example 
the presence of an equilibrium temperature gradient, a more general equation of state or diffuse interfaces, which 
might be important in certain applications^*^ However this study offers a systematic approach for examining the 
effects of compressibility which could represent a starting point for analyzing different or more complex configurations. 
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APPENDIX A: EQUATIONS FOR THE VISCOUS CASE 

Following the usual procedure for the incompressible, constant density case, the variables are nondimensionalized 

f 2 \ -1/3 / \ -1/3 

using 1/no = if— J as time scale and 1/fco = ( -jr\ as lengthscale. The compressibility parameter M is 

defined by M — ■ For simplicity the kinematic viscous coefficient is considered continuous over the interface, 

so that fJ.\/pi = /J.2/V2, with \i\ and /12 constant on each side of the interface. The value of the kinematic viscous 
coefficient at the interface is denoted by Voq. 

The scaled equations for u\ and A on each side of the interface can be written as 

A 4 D i u 1 + A 3 D 3 u 1 +A 2 D 2 u 1 +A 1 Du 1 +A Q u 1 =Q (Al) 
PiA = B 3 D 3 Ul + B 2 D 2 ui + BiDui + B (A2) 

where the coefficients (with the index m denoting the side of the interface suppressed for simplicity) are given by 



10 



A 4 = B 3 x 2 

A 3 = {DB 3 + B 2 )M2-B 3 co 

A 2 = (DB 2 +B 1 )p 1 [3 2 + exp(aMx)0 2 -B 2 uj 

A 1 = {DB\ + B o )0i02 — — — B\u> 

n 

A = DB o 1 2 - (n+ k 2 exp(aMx))0 2 - B uj 
B 3 



B 2 



eMaMx) ( 4 
aM exp(aMx) 



Pi 



R - " 

1 " ~4 



2 

r^,2 



7 



27 - 3 ) Ih 



a 2 M 2 expfaMx) I 1 . 
n h-l-3/% 



+ /?2 ( 7 + ^3) ( X + 77 exp(aMa;)^ 



B 



anM 



3 I 7 - 1 - i/3 3 ) + (27 - l)aMt 2 cxp(aMj)) 



Mk 2 2 



an* 



(A3) 
(A4) 
(A5) 

(A6) 

(A7) 

(A8) 

(A9) 



(A10) 
(All) 



DB 3 = 



~fa 2 M 2 ncxp{2aMx) 



a 3 m 3 nexp{2aMx) ( 67 - 4 
wd 2 = ^5 ry — d — p 3 



3/3| 



DBo = 



aMn 
a 2 M 2 n 



a 3 m 3 exp(2aMx) ^3 7 - 2 - l/3 3 ) + /? 3 2 /? 2 Q - 7^^) + /3f 
a 2 M 2 exp(2aM:r) 



/3 2 3 



-(7-ir (7-1- 3/33 



-(27 - \)aMk 2 exp(2aMx) 7 - -/?; 



3"°; ' aMn 2i2 



a 2 M 2 exp(aMx) 



n 



(7 - 1) ( 7 - 1 - 3^3 



Mfc 2 4fc 2 , A 
- /3 2 (l+7^ + ^exp(aM,)j 



02 = 7+3^3 

/3 3 = aMnexp(aMi) 
The equation for DA can be written as 



(A12) 
(A13) 
(A14) 

(A15) 

(A16) 

(A17) 
(A18) 



DA = 



aMn 
~02~ 



7 — 1 1 

A - exp(aMx)D 2 u 1 + -Du x + in + k 2 exp(aMx))ui 

n n 



(A19) 



while the equation for Z? 2 A is 
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D 2 A 



aMn 



aM(j - 1) 



P2 L n(3 2 
1 (27 - 1)M exp(aM:r) 



7 - 1 - -03 ) A - exp(aMa;)L'' :i wi 



'aM 
/ aMn 



1 



7 - 1 - -/3 3 ) + n + k exp(aMi) ) Dui 



\ 02 



7-1 



(27 - \)aMk 2 exp(aMi) ^ 
n/3 2 



Ul 



(A20) 



For 7 



00 or p Q 



00 (M -» 0), equations IA2I IA19I and IA201 yield A = 0, DA = 0, and D 2 A = 0, so the 



incompressible case is recovered. In the case 7 = 00 the equation for u\ simplifies to 



£> 4 ui - (nexp(-aMi) + 2k 2 )D 2 m + naM exp(-aMx)Dui 

nexp(-aMx) + exp(-aMx) + k 2 ) k 2 u x = (A21) 



If, furthermore, M — * in equation IA21I then the well known equation for uniform density incompressible fluid 
derived in Refi2 is obtained. The same equation can be obtained by letting p^ — > 00 (M — > 0) directly in equation 
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FIGURE CAPTIONS 

Fig. 1. Nondimensional rate of growth as function of jj = g (pf^° P2 ) f° r different values of 71 and 72. 

Fig. 2. Finite size effect on the nondimensional rate of growth, (a) No symbols curves represent the compressible 
case with 71 = 72 = 1.4, open symbols the corresponding constant density incompressible case and closed symbols 
the corresponding variable density incompressible case, (b) L\ = 0.1, L 2 = 0.05. All cases have At=0.5. 

Fig. 3. Compressibility influence on the nondimensional rate of growth in the presence of surface tension, for 
-, T f — r = 0.78 m 2 . All compressible cases have g ( pi+P2 ) = 4 mT 1 except the open symbol case for which g ( pi+P2 ) = 

g(Pl+P2) * Poo r Poo 

0.4 mT 1 . 

Fig. 4. Growth rate dependency on the the wave number for viscous fluids, (a) The compressible and incompress- 
ible variable density cases have M = 0.1. (b) At=0.1, compressible and variable density incompressible cases showed 
with thick lines correspond to M = 0.1 and with thin lines to M = 0.01. 
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